Reynolds numbers of the large-scale flow in turbulent Rayleigh-Benard convection 
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We measured Reynolds numbers Re of turbulent Rayleigh-Benard convection over the Rayleigh- 
number range 2 x 10^ ^ R ^ 10^^ and Prandtl-number range 3.3 < a < 29 for cylindrical samples of 
aspect ratio T = 1. For i? ^ 7?c ~ 3 x 10^ we found Re - R'^"" with /J^// ^ 0.46 < 1/2. Here both 
the a- and 7i-dependences are quantitatively consistent with the Grossmann-Lohse (GL) prediction. 
For R> Rc we found Re = 0.106 a'^^'^R^^^, which differs from the GL prediction. The relatively 
sharp transition at Rc to the large- J? regime suggests a qualitative and sudden change that renders 
the GL prediction inapplicable. 
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Understanding turbulent Rayleigh-Benard convection 
(RBC) in a fluid heated from below 1] remains one of 
the challenging problems in nonlinear physics. It is well 
established that a major component of the dynamics of 
this system is a large-scale circulation (LSC) lit. For 
cylindrical samples of aspect ratio T = D / L :^ 1 {D is the 
diameter and L the height) the LSC consists of a single 
convection roll, with both down-flow and up-flow near the 
side wall but at azimuthal locations 9 that differ by tt. An 
additional important component of the dynamics is the 
generation of localized volumes of relatively hot or cold 
fluid, known as "plumes" , at a bottom or top thermal 
boundary layer. The hot (cold) plumes are carried by 
the LSC from the bottom (top) to the top (bottom) of 
the sample and by virtue of their buoyancy contribute to 
the maintenance of the LSC. The LSC plays an important 
role in many natural phenomena, including atmospheric 
and oceanic convection, and convection in the outer core 
of the Earth where it is believed to be responsible for the 
generation of the magnetic fleld. In this Letter we report 
measurements of the speed of the LSC that agree well 
with a theoretical prediction by Grossmann and Lohse 
for relatively small applied temperature differences 
AT, but depart from this prediction rather suddenly as 
AT is increased further. Our results illustrate clearly 
that a quantitative understanding of this system is still 
restricted to limited parameter ranges. 

The LSC can be characterized by a turnover time T 
and an associated Reynolds number 4] 



Re = (2T/r) X (L/iy) 



(1) 



(v is the kinematic viscosity). A central prediction of var- 
ious theoretical models P, S 0, ^ is the dependence 
of Re{R,(j) on the Rayleigh number 



R = agATL^/Ki^ 



and on the Prandtl number 



a = V I K, 



(2) 



(3) 



A recent prediction by Grossmann and Lohse (GL) 
based on the decomposition of the kinetic and the ther- 
mal dissipation into boundary-layer and bulk contribu- 
tions, has been in remarkably good agreement with ex- 
perimental results for Re{R,a) However, the param- 
eter range covered by the measurements was relatively 
small. 

We report new measurements of Re{R, a) over a wider 
range, for R up to 10^^ and 3.3 < cr < 29. For modest 
R, say R < 2 X 10^, we again find very good agreement 
with the predictions of GL. However, for larger R the 
measurements reveal a relatively sudden transition to a 
new state of the system, with a Reynolds number that is 
described well by 



R, = 0.106 a-^/^ R^'^ 



(4) 



This result differs both in the cr dependence and in the 
R dependence from the GL prediction. We interpret our 
results to indicate the existence of a new LSC state. It 
is unclear at present whether the difference between this 
state and the one at smaller R will be found in the ge- 
ometry of the flow, in the nature of the viscous boundary 
layers that interact with it, or in the nature and frequency 
of plume shedding by the thermal boundary layers adja- 
cent to the top and bottom plates. But whatever its 
nature, this state does not conform to the consequences 
of the assumptions made in the GL model. 

Another important aspect of the predictions is the de- 
pendence of the Nusselt number (the dimensionless effec- 
tive thermal conductivity) 



M=QL/XAT 



(5) 



(a is the isobaric thermal expansion coefficient, n the 
thermal diffusivity, and g the acceleration of gravity). 



on R and a (here Q is the heat-current density and A 
the thermal conductivity) . The GL model 0, 0| provides 
a good fit also to data for M at modest R, say up to 
Here we briefly mention as well 
measurements of Af for larger R jl2] that depart signifl- 
cantly from the GL prediction as R approaches 10^^. 

Measurements of Re were made for three cylindrical 
samples with F ~ 1. Two of them, known as the medium 
and large sample, 0| had L — 24.76 and 49.69 cm re- 
spectively. The third was similar to the small sample of 
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FIG. 1: The cross-correlation function (squares) between two 
thermometers (1 and 5) mounted on opposite sides of the side 
wall, and the auto-correlation function (circles) of a single 
thermometer (1), on a logarithmic scale as a function of the 
delay time r for the large sample and R — 7.5 x 10^". Solid 
lines: fits of Eqs.E|or|S]to the data. The lengths of the lines 
indicate the range of the data used in the fits. 



Rcf. 13], but had L = 9.52. As evident from Eq. 13 a 
given accessible range of AT will provide data over dif- 
ferent ranges of R for the different L values. For the 
small sample we used 2-propanol with cr = 28.9 as the 
fluid and measured the frequency / of oscillations of the 
direction of motion of plumes across the bottom plate 
to obtain i?e = 21"^//^ |3. With the medium and 
large sample we used water, mostly at mean temper- 
atures Tm = 55.00, 40.00, and 29.00°C corresponding 
to (T = 3.32, 4.38, 5.55 and v = 5.11 x IQ-"^, 6.69 x 
10"'', 8.25 X 10"'' m^/sec respectively. The top and 
bottom plates were made of copper. A plexiglas side 
wall had a thickness of 0.32 (0.63) cm for the medium 
(large) sample. At the horizontal mid-plane eight ther- 
mistors, equally spaced around the circumference and la- 
beled i = 0, ...,7, were imbedded in small holes drilled 
horizontally into but not penetrating the side wall. The 
thermistors were able to sense the adjacent fluid temper- 
ature without interfering with delicate fluid-flow struc- 
tures. When a given thermistor (say i = 0) sensed a 
relatively high temperature Ti due to warm upflow of 
the LSC, then the one located on the opposite side (say 
at i = 4) would sense a relatively low temperature due 
to the relatively cold downflow. 

When a warm (cold) plume passed a given side-wall 
thermistor, the indicated temperature was relatively high 
(low). It had been shown before by comparison of 
temperature sensors actually imbedded in the fluid and 
laser-doppler velocimetry, that this thermal signature can 
be used to determine the speed, and thus the Reynolds 
number, of the LSC and that it yields the same result 
as actual velocity measurements. Indeed, where there is 
overlap, our results for Re are in satisfactory agreement 
with measurements Q based on velocimetry. 

From time series of the eight temperatures Ti(t) taken 
at intervals of a few seconds and covering at least one 



and in some cases more than ten days at each of many 
values of R we determined the auto-correlation functions 
(AC) Cij{T), i = j, and the cross-correlation functions 
(CC) Cij{T), i = 0, . . . ,3, j — i + i corresponding to 
signals at azimuthal positions displaced around the circle 
by TT. They are given by 



C.Ar) = mit)-mt))t] X [T,{t + r)-{T,it))t])t . (6) 

We show an example of AC (circles) and of CC (squares) 
in Fig.n 

One sees that the AC have a peak centered at the ori- 
gin. It can be represented well by a Gaussian function. 
The peak width indicates that the plume signal is corre- 
lated over a significant time interval. A second smaller 
Gaussian peak is observed at a later time ij'^ that we 
identify with one turn-over time T of the LSC. The exis- 
tence of this peak indicates that the plume signal retains 
some coherence while the LSC undergoes a complete ro- 
tation 15]. A further very faint peak is found at 2T, 
but is not used in our analysis. These observations are 
consistent with previous experiments I1E1E3- This 
structure is superimposed onto a broad background that 
decays roughly exponentially on a time scale of O(IOT). 
We believe that the background decay is caused by a slow 
meandering of the azimuthal orientation of the LSC. 

The CC are consistent with the AC. Here too there is 
a broad, roughly exponential, background. There is no 
peak at the origin, and the first peak, of Gaussian shape, 
occurs at a time delay tf^ = T/2 associated with half a 
rotation of the LSC. A further peak is observed at 3T/2, 
corresponding to 1.5 full rotations. 

Based on the above, we fitted the equation 



Ci,iiT) = boexp ( -— ) +biexp 



-ac 
'l 



b2exp 



(7) 



to the data for the AC, and the equation 



-biexp 



(8) 



to those for the CC. Examples of the fits are shown in 
Fig. n as solid lines. One sees that the fits are excellent. 
Substituting T = and T ^ 2tf into Eq.[Tl we have 
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FIG. 2: (a): Reynolds numbers Rf^ for a = 4.38 as a function 
of the Rayleigh number _R on logarithmic scales, (b): Re/R^^^ 
on a linear scale, as a function of R on logarithmic scales. 
From top to bottom, the four data sets are for a = 3.32, 
4.38, 5.55, and 28.93. Solid (open) circles: RT {RT), large 
sample. Solid squares: i?^*^, medium sample. Solid diamonds: 
-Re for 2-propanol at 40°C. Open triangles: Ref. cr = 
5.4. Dashed lines (from top to bottom): GL predictions for 
a = 3.32, 4.38, 5.55, and 28.9. O Horizontal solid lines: 
R^/R^'^ = 0.0433, 0.0347, and 0.0293. 



as two experimental estimates of Re- For each R we 
computed the average value i?^'^ of the eight CC Cij and 
Cj^i with i — 0, . . . , 3 and j = i + 4. The results are 
shown in Fig. |21 as solid squares (medium sample) and 
solid circles (large sample). Averages of the eight AC 
Ci^i at each R for the large sample are shown as open 
circles. There is excellent agreement between the AC and 
the CC. Also shown, as solid diamonds, are results for the 
small sample deduced from the oscillation of the direction 
of plume motion across the bottom plate Q. These 
data are for 2-propanol with a = 28.9. For comparison, 
the results of Qiu and Tong 0, 0| based on velocity 
measurements for a = 5.4 are shown as open triangles. 
Our data for a = 5.55 are in quite good agreement with 
them. 

The dashed lines in Fi g. El are, from top to bottom, 
the predictions of GL Ulla for a = 3.25, 4.38, 5.55 
and 28.9. For i? < 3 x 10^ they pass very well through 
the data. We regard this agreement of the prediction 
with our measurements as a major success of the model. 
However, for larger R the data quite suddenly depart 




FIG. 3: Effective exponents /9e// (a) and Qe// (b) of Eg. 1111 
as a function of R. Dotted lines: asymptotic values predicted 
by GL for R oo. Dashed lines: effective values as a func- 
tion of R predicted by GL for a — 4.38. Solid circles in (a): 
experimental values for /3e// obtained by fitting a powerlaw 
to the data for Re{R) at ct = 4.38, using a sliding window 0.8 
decades wide. Dash-dotted line in (a): /3e// ~ 1/2. Solid line 
in (b): approximate location of experimental results. 
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FIG. 4: The reduced Reynolds number R^a^^'^ R'^^'^ as a 
function of the Prandtl number a on logarithmic scales. 
Dashed lines (from top to bottom): GL predictions for R = 
" "ill- Solid line: Eq.|3| Dots: all data for 



10^ 10^°, and lO" 
3 X 10^ < R < 10" 



from the prediction and scatter randomly about the hor- 
izontal solid lines. These results indicate that there is a 
sudden change of the exponent /3e/ / of the power law 



Re{R, o) = RqG 



(11) 



as R exceeds i?c — 2 x 10^, from a value less than 1/2 
to 1/2 within experimental resolution. The GL model 
can not reproduce this behavior, and we conclude that 
a new large- i? state is entered that does not conform to 
the assumptions made in the model. 

The inconsistency between the prediction and the data 
can bee seen more clearly by considering the effective 
exponents /3e// and cHe// defined by Ea. llll and shown in 
Fig. El The experimental values of (3ef f in Fig. |3Ji were 
obtained by fitting powerlaws to the data for a — 4.38, 
using a sHding window 0.8 decades wide. One sees that, 
within experimental uncertainty, the value 1/2 is reached 
at i? ~ 7 X 10^. As can be seen from Fig. [3 Peff ^1/2 
actually is reached earlier, near i? ~ 3 X 10^; the results 
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agree well with the prediction by Grossmann and Lohse 
0; but for larger R we find that Re = 0.106a^/^R^/^ , in 
disagreement with the GL prediction. 

Wc thank Xin-Liang Qiu for providing us with the nu- 
merical data corresponding to Fig. 12 of Ref. This 
work was supported by the US Department of Energy 
through Grant DE-FG02-03ER46080. 



FIG. 5: The reduced Nussel number Nao/B}'^^ as a function 
of -R for cr = 4.38 Dashed line: GL prediction. 

in Fig. |3t represent an average over a finite range of R 
because a finite window width had to be employed in the 
analysis. The GL model predicts the value /3 = 4/9 ~ 
0.444 when R becomes large enough so that a pure power 
law prevails (dotted line). The predicted effective values 
fieff at finite R (dashed line) are already very close to 
this value in the experimental range of R. It is hard 
to see how the prediction could be changed by adjusting 
parameters in the model so as to yield j3ef / = 1/2 for R > 
2 X 10^ without changing the seemingly firm prediction 
(3f.ff ^4/9 for sufficiently large R. 

Figure shows ctefj defined by Eq. ^2 Here the 
GL prediction yields a = 2/3 as i? becomes large (dot- 
ted line). The dashed line gives the predicted ae// as 
a function of R for finite R. One sees that cte// is al- 
ready quite close to a in the experimental range of R. 
In the range where the experimental Re/R^^^ is constant 
(i.e. 2 X 10^ < R < 10^^) the data are consistent with 
aeff = 3/4 as shown by the solid line, but not with 
aeff — 2/3. To explore this point further, we show in 
Fig.0]i?eCr'^/'*/i?^/^ as a function of a. Here all our data 
for i? > 3 X 10^ arc plotted. We note that, at each 
cr, all data collapse into a narrow range consistent with 
the scatter of the measurements. The horizontal solid 
line, which corresponds to Eq.^and thus to aeff = 3/4, 
passes through the points at each cr within the scatter, 
although it is a bit low for the smallest a. The solid 
lines are the GL predictions for, from top to bottom, 
R = 10^ 10^°, and lO". 

It is interesting to note that a similar inconsistency was 
found also between the GL prediction and the Nusselt 
number This is illustrated in Fig. |31 where we show 
the reduced Nusselt number J\f^/R^^^ as a function of 
R. There are deviations from the prediction (dashed 
line) for R > 10^°, which is somewhat higher than the 
value of Rc for the Reynolds number. 

In this Letter we presented new measurements of the 
Reynolds number Re of the large-scale circulation in tur- 
bulent Rayleigh-Benard convection for an aspect-ratio- 
one cylindrical sample over the Rayleigh-number range 
2 X 10^ < i? < 10" and the Prandtl-number range 3.3 < 
cr < 29. For i? ^ 3 X 10^, where Re oc R°-^^, our data 
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